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We discuss a combination of unitary and anti-unitary symmetry of quantum Liouvillian dynam- 
ics, in the context of open quantum systems, which implies a D2 symmetry of the complex Liou- 
villian spectrum. For sufficiently weak system-bath coupling it implies a uniform decay rate for all 
coherences, i.e., off-diagonal elements of the system's density matrix taken in the eigenbasis of the 
Hamiltonian. As an example we discuss symmetrically boundary driven open XX Z spin 1/2 chains. 
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Introduction. - The so-called PT symmetry, a product 
of a unitary and an anti-unitary transformation both of 
which square to identity, of non-Hermitian Hamiltoni- 
ans has been introduced by Carl Bender [T] to study 
the spectral theory of Schrodinger-like operators and 
non-standard formulations of quantum mechanics. The 
generic possibilities ^-[5] of having non-Hermitian opera- 
tors with purely real spectra have recently found impres- 
sive experimental applications in non-linear optics [71 |S] 
and even in LRC electric circuits ^ where PT symmetry 
is achieved by a subtle combination of active elements, 
with the gain and the loss distributed in a symmetric 
way. 

In this Letter we show that the concept of PT symme- 
try can be introduced, in contrast to Ref. Ij, in the con- 
text of standard, orthodox quantum mechanics when one 
considers open, dissipative systems. We discuss a gen- 
eral situation in which a Liouvillian superoperator pos- 
sesses a combination of unitary and anti-unitary master- 
symmetries (transformations in the linear space of gener- 
ators of hermiticity preserving dynamical semi-groups). 
We show that, as a consequence of master symmetry, 
the spectrum of decay rates should have a dihedral (-D2) 
symmetry in the complex plane. Furthermore, for suffi- 
ciently weak system-bath coupling the spectrum has the 
shape of a cross, with one leg on the real axis, correspond- 
ing to dynamics of populations (diagonal operators) in 
the energy eigenbasis, and the other leg parallel to the 
imaginary axis, corresponding to the decay of coherences 
(off-diagonal operators), remarkably, all with the same 
asymptotic rate. This phenomenon offers a fundamen- 
tally new way of controlling decoherence as one deals 
with a single damping factor. Therefore, our general re- 
sult should be of interest for a variety of fields which use 
the methods of open quantum systems [10] , ranging from 
non-equilibrium statistical mechanics, quantum optics, 
quantum information, and quantum measurement the- 
ory, to condensed matter, high-energy theory and quan- 
tum cosmology. 

As an example, we demonstrate how our symmetry is 
realized in the symmetrically boundary driven XX Z spin 
1 /2 chain [TTI - [l4] which is described in terms of the stan- 
dard Lindblad equation (with Hermitian Hamiltonian) 



and the canonical formalism of Markovian open quan- 
tum systems. We stress that our considerations assume 
only local-and-homogeneous-in-time open quantum sys- 
tem's theory, without any need for non-Hermitian central 
system's Hamiltonians. Yet, the results we find formally 
faithfully generalize the mathematical framework of PT- 
symmetric quantum mechanics [J to a Liouvillian set- 
ting. 

Quantum dynamics and Liouville space formalism. - 
Let us consider a system defined on a finite Hilbert space 
H oi N states, with a canonical orthonormal (ON) ba- 
sis \j},j = 1 ■ • ■ N. Let B{T-L) denote the vector space of 
linear operators over H. N'^ basis states of B{T-L) shall 
be denoted by Ej^k — j,k = 1...N. Introduc- 

ing the Hilbert-Schmidt inner product (see e.g. [T5] ) 
(cr, p) :— tr(crV), BCH) becomes a Hilbert space, with 
{Ej^k} being its ON basis. As arbitrary physical states 
are elements of B{H) in the sense of density operators, 
the generator C of quantum Liouvillian dynamics 

^^p{t)^Cp{t) (1) 

can be considered as an element of B{B{H)) and can 
be, say in a basis {Ej ^}, represented by iV^ x ma- 
trix. For general Markovian open quantum systems, the 
Liouvillian C can be always cast to the Lindblad form 

[inidSKTi] 

C = -i(ad H) + -f-D, (ad H)p := [F, p] , 
'Dp ■■= ^ 2L„,pLl^ - Ll,L„,p - pLl^L,^, (2) 

rn 

separating the unitary part generated by a Hermi- 
tian Hamiltonian H e B{H), and the dissipator V e 
B{B{H)), where {L,„;m = 1 . . . M} C B{U) is a set of 
M < N"^ — 1 Lindblad operators, which together with the 
system-bath coupling strength 7 > contain all informa- 
tion that is left about the reservoirs (i.e., environment 
degrees of freedom) and coupling to them. One can al- 
ways adjust 7 such as to fix the trace of the dissipator 
Tri* = T,j,kiEj,k,T>Ej,k) = - Tri = -N^. We note 
that the abstract part of our discussion will not require 
a particular Lindblad form ^ of the Liouvillian C, but 
we require only the existence of the following concepts. 
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Liouvillian -transformation.- Let V €E B{B{H)) be 
some unitary transformation, ("Per, Vp) = {a, p) , Vu, p £ 
B{'H), squaring to identity = 1. Noting that V^^ ~ 
V, we define P-transformation of the Liouvihian as 



PC := VCV. 



(3) 



Liouvillian T -transformation.- The Hilbert-Schmidt 
inner product completely defines the Hermitian adjoint of 
the LiouvUean C\ namely, {a,C^p) = {Ca, p) — {p,Ca-), 
y<J,p e B{'H). We define the T-transformation as 



(4) 



Note that the maps P and T can be considered as ele- 
ments of B{B{B{7i))). By introducing the supcr-Hilbert- 
Schmidt inner product in {{X,y)) TriX^^y), 

it can be verified straightforwardly that the map P is uni- 
tary, while the map T is anti-unitary, namely, 

{{vx,vy)) = {ix,y)), (5) 
{iTA:,Ty)) = {iy,x)), yxje B{B{n)). (6) 

¥T- symmetric Liouvillian.- Let C denote the traceless 
part of the Liouvillian, defined as 



£ = £'-71, with 7:=-Tr£/Tri. 



(7) 



We shall define quantum Liouvillian dynamics as PT- 
symmetric if the following key identity holds 



¥TC' 



(8) 



or equivalently, {C'Y = —VC'V. This PT symmetry can 
be immediately applied to invert the Liouvillian propa- 
gator of generally dissipative dynamics hl{t) :— exp(<>C) 



U{-t) = e'^'V[U{t)\^V. 



(9) 



Now we are in position to show the following result: 
Theorem. The spectrum of WT- symmetric Liouvillian 
has a dihedral group (D2) symmetry in the complex 
plane, with the lines of symmetry — —7 + iM and 

= R- More specifically, writing the Liouvillian spectral 
decompositions 



(10) 



with right and left eigenvectors which can be chosen bi- 
orthonormal, {ua,vp) = Sa^p, we .state that for each 
eigenvalue of C , = -f- 7, there exist the eigenval- 
ues X'p = —X'a (image across £v), and X'^ = X'a (image 
across ih) with the right and left eigenvectors related via 



(11) 

(12) 



In the case of a degenerate eigenvalue Aq,, the eigenvec- 
tors upf^ and vp/r^ are to be understood as appropriate 
members of the right and left eigenspaces, respectively. 



Proof. The fv-reflection spectral symmetry with (11) 
is a direct consequence of PT symmetry ([s] ), after ap- 
plying it to the left-hand-sides of Eqs. (10), and mul- 
tiplying the resulting equation by V. The £h-reflection 
symmetry with (|12|), on the other hand, is an immedi- 



ate consequence of hermiticity preservation of Liouvillian 
quantum dynamics, C{p^) = {Cp)\ which clearly holds 
for the Lindbladian ([2]), but also for any other (possibly 
non-Markovian) meaningful quantum dynamics |17j . □ 

Remarks. - One of the most interesting objects of open 
quantum dynamics is the steady state, p{t — 00) = ui 
corresponding to eigenvalue Ai = 0, which always exists, 
due to trace preservation in %, which can be expressed 
as C^vi = with vi = 1 = X^jLi |j)OI- This means that 



Eq. (11) yields already a non-trivial result, namely that 
the fastest decaying mode A^v^ = — 27 is Mjvs = ■p(l). 

Even more remarkable observation is that for suffi- 
ciently small coupling 7, say below some critical value 
7 < 7pT) the spectrum of decay modes lies strictly on the 
cross, {Aq,} C fvUfh- First, we show that if an eigenvalue 
Aq(7) lies on the real line £h, it remains on as long 
as the eigenvalue is isolated. First order non-degenerate 
perturbation theory tells us that dAQ/d7 = (wQ,I?Ma), 
whence hermiticity conservation VjugY = 2?(ujj,), and 



Va following from (12 1 since Aq = Aq, im 



plies dAQ/d7 = dAQ/d7 g M. Second, we show similarly 
that for an isolated eigenvalue Aq(7) initially on i^, i.e., 
Aq G iE, we have dAQ/d7 g iM as a simple consequence 
of the PT symmetry of the traceless part of the dissipator 
X>' = P + i, VTV' = -V'. Namely^ from = -A^ and 
(11) follows that Mq = Vva, so dA'a/d7 = (2?'uq,Uq) = 
-{uo^VV'Vva) = -(wq,P'mq) = -dA'Q/d7. The spec- 
trum can then leave the cross £v U ^h, when at some 
7 = 7PT a pair of eigenvalues collides and shoots off 
into the complex plane. By observing just the motion 
of the spectral points on the vertical leg (while the 
argument should be quite similar for £h) the critical cou- 
pling strength where the PT symmetry of the spectrum 
is spontaneously broken can be estimated heuristically as 



7PT 



\V'\ 



-1 J-2 



(13) 



Here ||2?'|| is the operator norm of the dissipator which 
estimates the maximal velocities |dAQ/d7| and d denotes 
a typical density of states of H , d^ giving a typical density 
of energy differences ej — t^- 

What remains to be shown to prove this picture is that 
initially, for (infinitesimally) small 7, the eigenvalues Aq 
indeed start on M U iM. This is again easy to demon- 
strate, working in the eigenbasis of the Hamiltonian H , 
H\i^j) = J = 1 . . . iV. For 7 = 0, A^, = i(e, - e^) 

where a now labels all N"^ pairs of j, fc. Let us assume 
that the energy spectrum £_,- is non-degenerate, i.e., all 
tj are different. Then we have exactly N diagonal eigen- 
operators dj — = "q = Wq for which Aq = 0. 

In order to understand their motion as we switch on 
7, we have to solve first-order-degenerate-perturbation 



3 



problem; i.e., we have to diagonalize an N x N matrix 
Vj-.fc := (djjVdk), Vua = S.ada, where each eigenvalue 
determines the motion of some A^, dA„/d7|^=o = Cqi ^^'^ 
the corresponding eigenvector determines the hybridiza- 
tion, Ua = Va = ^ciajdj. The key observation now is 
that the matrix V is real and symmetric [18j 



Fj-fc - idj,v'dk) = {^p,\v'i\i>k){4'k\M) 

+ Sj,k = Vj^k = Vk^j, (14) 

proving that these N Liouvillian eigenvalues Xail) 
main on the real line. Asymptotically, for small 7 the 
states ^0,(7) are diagonal and correspond to populations 
in the energy eigenbasis. The other N'^ — N eigenvalues 
A^(7), for ej ^ efe, move on the imaginary line, as shown 
in the previous paragraph, provided that they are iso- 
lated initially for 7 = 0, i.e., provided that all the energy 
spacings ej — Cfc are different (non-degenerate). These 
eigenvalues, asymptotically for small 7, correspond to off- 
diagonal operators - coherences - in the energy eigenba- 
sis, and all decay with exactly the same rate RcAq = —7. 
This is remarkable and should have experimentally ob- 
servable consequences, e.g. one should be able to control 
the decoherence in such a system by handling a single 
damping factor e~'^*. 

The above scenario which predicts the full Liouvil- 
lian spectrum of decay modes to belong to the cross 
£v U £h for some non-empty coupling strength interval 
< 7 < 7pt is strictly justified only if the two condi- 
tions are met, namely that both, the energy spectrum 
{ej} and the energy difference spectrum {ej — ek',j 7^ k} 
are non-degenerate. In the non-generic case when we 
have a degeneracy in either of the two, it can happen 
(in the absence of additional selection rules) that already 
infinitesimal dissipation 7 moves the corresponding Li- 
ouvillian eigenvalues out into the complex plane, so the 
spontaneous PT symmetry breaking of the Liouvillian 
spectrum may then occur already for a vanishing per- 
turbation 7pt — 0. 

Example.- We close by demonstrating our construc- 
tions in an interesting example, namely we consider an 
open XX Z chain of n spins 1/2 with the Hamiltonian 

n-l 

^ = E (2^>7+i + 2^7 <+i + (15) 

where = ^(^j ^ = 1 ■ • • ?t- are Pauli oper- 

ators on a product space H ~ (C^)*^", with symmetric 
Lindblad driving acting on the edges of the chain only 

il,2 = ^VT±~/i erf , ^3,4 = ^Vl TMCT,f , (16) 

with M = 4, where G [—1,1] is a driving parameter 
determining the magnetization bias between the left and 
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FIG. 1. The Liouvillian spectrum {A^} for an open XXZ 
chain with A = 1/2, n = 4, maximum driving ^ = 1, and 
7 — 0.02 < 7pt (top), 7 = 0.2 > 7pt (middle), and 7 = 2 > 
7pt (bottom) . Note the dihedral symmetry of the spectrum 
with the horizontal and vertical symmetry lines, £h = K and 
£v ~ —7 + iK, respectively. Because the XXZ chain has a 
global conservation law — '^jC7j, we consider only the 
most relevant sector with zero total magnetization — 0. 



the right bath. We have here N — This model has 
been intensively studied recently [TTUni [TU] , and it has 
been shown to admit exact solutions [14) in the limiting 
cases of small 7 or /i = 1, and can exhibit diffusive spin 
transport (for |A| > 1) in the linear response regime |13j . 
of small fi and non-small 7 ~ 1. In order to disclose ex- 
plicitly the PT symmetry of such a symmetrically bound- 
ary driven XXZ chain, it is instructive to identify the 
operator space B{H) with the tensor product H(E)'H, via 
the isomorphism 



|V')(</'|o|i^)®5|<^) 



(17) 



where S :~ YYj=i '^^ a global spin-flip operation in (t| 
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eigenbasis. Then the Liouvilhan is represented as 



7(1 + m) 
2 

7(1 -m) 



{(7^ ®CJ^ + CJ,, ® cr+) 

{a^ + a+®(7~)--fl®l. (18) 



Let us define the parity transformation V, such that it 
corresponds to the foUowing operator 



chain [T^ [20] is PT— symmetric as well (as can be eas- 
ily seen in spin-ladder formulation), but applications to 
driven open bosonic cold atom systems should also be 
possible. 
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R 



(19) 



where i? is a reflection permutation which reverses the 
order of sites j < — > 
mon eigenbasis of a^. 

X]mi...m„G{ + l,-l} 1"^" ■ 

note U = i?n, '^? 



n + 1 — j, i.e., in the com- 
j = 1 . . .n, it reads R = 
. mi)(r7ii . . . m„|. In notation of 
, W = R. It takes a straightfor- 
ward calculation to show that indeed 7^^ = 1 O 1 1, 
and (£ + 71)1' = -ViC + -fi)'P. In fact, the PT sym- 
metry |8]) holds for each of the three rows of the expres- 
sion ( |18[ ) separately. In Fig. [T] we show three different 
Liouvillian spectra for different values of the coupling 
constant, before and after the transition 7 = 7pt. Nu- 
merical experiments indicate that the critical value in 
the leading order decays exponentially with the chain 
length, 7pt oc oc 4~", as estimated in (13 1. We 



note that even for 7 3> 7pt, a substantial fraction of 
spectral points remain on the line £v = ^7 + ilR, hence a 
significant spectral weight for uniform relaxation e~^* is 
expected for typical (offdiagonal) observables. Remark- 
ably, an important operator in the transport theory, the 

spin current operator J — iJ2]j=ii'^j''^J+i~'^j''^^+i)^ 
vanishing diagonal matrix elements in the energy eigen- 
basis, so its expectation value relaxes with a uniform rate 
I tr [(pit) - p(oo)) J]| ~ e-^*, for 7 < 7pt. 

Discussion. ~ We outlined a general framework for anal- 
ysis of a combined unitary (P) and anti-unitary (T) 
master-symmetry of the most general types of quantum 
master equations which are local in time. We stress that 
our analysis remains strictly in the framework of canoni- 
cal quantum mechanics, so we need no active elements or 
non-Hermitian system's Hamiltonians for our construc- 
tions. PT symmetry of dissipative Liouvillian dynamics 
can thus occur only with respect to a shift parallel to 
the imaginary line which represents an average damping 
rate. In the asymptotic regime of weak system-bath cou- 
pling, the Lioivillean spectrum can be strictly separated 
with the coherences - the off-diagonal matrix elements 
of the state in the energy eigenbasis - decaying with 
a strictly uniform rate. We discussed a simple explicit 
example of PT-symmetric Liouvillian dynamics in open 
XX Z spin chains, but other interesting and experimen- 
tally accessible realizations are possible. For example, 
the recently studied symmetrically driven fermi Hubbard 
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